Statistical Determination of Bulk Flow Motions 
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We present here a new parameterization for the bulk motions of galaxies and clusters (in the linear 
regime) that can be measured statistically from the shape and amplitude of the two-dimensional 
two-point correlation function. We further propose the one-dimensional velocity dispersion (v p ) of 
the bulk flow as a complementary measure of redshift-space distortions, which is model-independent 
and not dependent on the normalisation method. As a demonstration, we have applied our new 
methodology to the C4 cluster catalogue constructed from Data Release Three (DR3) of the Sloan 
Digital Sky Survey. We find v p = 270 +433 km/s (also consistent with v p = 0) for this cluster sample 
(at z = 0.1), which is in agreement with that predicted for a WMAP5-normalised ACDM model 
(i.e., t> p (ACDM) = 203 km/s). This measurement does not lend support to recent claims of excessive 
bulk motions (~ 1000 km/s) which appear in conflict with ACDM, although our large statistical 
error cannot rule them out. From the measured coherent evolution of v p , we develop a technique to 
re-construct the perturbed potential, as well as estimating the unbiased matter density fluctuations 
and scale-independent bias. 

PACS numbers: draft 



I. INTRODUCTION 

A decade ago, astronomers discovered the expansion 
of the Universe was accelerating via the cosmological 
dimming of distant supernovae P, HJ. Since then, the 
combination of numerous, and diverse, experiments has 
helped to establish the Cosmological Constant (specifi- 
cally a ACDM model) as the leading candidate to explain 
this cosmic acceleration. However, with no theoretical 
motivation to explain the required low energy vacuum of 
the ACDM model, there is no reason to preclude alter- 
native models, especially those based upon the possible 
violation of fundamental physics which have yet to be 
proven on cosmological scales 0, H[ . 

In addition to using geometrical probes like Super- 
novae to constrain the cosmic acceleration, tests based 
on the formation of structures in the Universe also pro- 
vide a method for validating our cosmological models. 
In particular, we can investigate the consistency between 
the geometrical expansion history of the Universe and 
the evolution of local density inhomogencities to help re- 
veal a deeper understanding of the nature of the cosmic 
acceleration 

In general, there are three observables that can be used 
to quantify structure formation in the Universe, namely 
geometrical perturbations, energy-momentum fluctua- 
tions and peculiar velocities, all of which will be measured 
to high precision via future experiments like DES, LSST, 
JDEM and Euclid (see details of these experiments in the 
recent FoMSWG report [l(| ) ■ In more detail, such weak 
lensing experiments measure the integrated geometrical 
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effect on light as its trajectory is bent by the gravitational 
potential. Likewise, galaxies (and clusters of galaxies) 
measure the correlations amongst large-scale local inho- 
mogencities, while the observed distortions in these cor- 
relations (in redshift-space) can be used to extract infor- 
mation about peculiar velocities (nl - flij . In this paper, 
we explore the cosmological constraints on the physics of 
cosmic acceleration using peculiar velocities, as it is one 
of the key quantities required for a consistency test of 
General Relativity [II [lj]. 

Early observational studies of the peculiar velocity 
field, or "bulk flows" , have produced for many years dis- 
crepant results [13], primarily due to small sample sizes 
and the heterogeneous selection of galaxies. However, a 
recent re-analysis of these earlier surveys fl8j has now 
provided a consistent observational picture from these 
data and finds significant evidence for a larger than ex- 
pected bulk motion. This is consistent with new mea- 
surements of the bulk motion of clusters of gala xies us- 
ing a completely different methodology (l9l |20| . which 
leads to the intriguing situation that all these measure- 
ments appear to be significantly greater in amplitude, 
and scale, than expected in a concordance, WMAP5- 
normaliscd ACDM cosmological model. Such discrepan- 
cies with ACDM may give support to exotic cosmological 
models like modified gravity [2l| . 

Given the importance of these large-scale bulk flow 
measurements, we propose here an alternative methodol- 
ogy to help check these recent claims of anomalously high 
peculiar velocities which are inconsistent with the stan- 
dard ACDM cosmology. We start by outlining a statis- 
tical determination of bulk flow motions using redshift- 
space distortions in large-scale galaxy or cluster surveys. 
Such redshift-space distortions are easily seen in the two- 
dimensional correlation function (£ s (c, 7r)), which is the 
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decomposition of the correlation function into two vec- 
tors; one parallel (tt) to the line-of-sight and the other 
perpendicular (a) to the line-of-sight. On small scales, 
any incoherent velocities of galaxies within a single dark 
matter halo (or cluster) will just add to the cosmological 
Hubble flow thus causing the famous "Fingers-of-God" 
(FoG) effect which stretches the 2-D correlation function 
preferentially in line-of-sight (tt) direction. These dis- 
tortions depend on the inner dynamics and structure of 
halos and therefore, any cosmological information is dif- 
ficult to distinguish from the halo properties. However, 
on large scales (outside individual dark matter halos), the 
peculiar velocities become coherent and follow the linear 
motion of the matter thus providing crucial information 
on the formation of large-scale structure [22j ■ 

In this paper, we compare predictions for £ s (tt, a) 
to observations based on the C4 cluster catalogue [23j 
from the Sloan Digital Sky Survey (SDSS) dj. Us- 
ing the Kaiser formulation [221 ]. a theoretical model for 
£ s (tt, a) is fit to the measured 2-D correlation function 
in configuration-space, with the £ s parameterised by a 
shape dependent part and a coherent evolution compo- 
nent. We also propose that the 1-D linear velocity disper- 
sion (v p ) is a interesting quantity to report when measur- 
ing redshift-spacc distortions, and complementary to tra- 
ditional quantities like /?, / or fag discussed recently [fll - 
[l3| , as it is independent of both bias and the normalisa- 
tion method. Therefore, the measured v p provides an 
unbiased tracer of the evolution of structure formation. 



II. STATISTICAL DETERMINATION OF 
PECULIAR VELOCITY 

The redshift-space two-point correlation function of 
mass tracers (£s(cr, tt)) is an anisotropic function [22j]. 
On small scales, it is elongated in the 7r-direction by the 
"Fingcrs-of-God" effect, while on large scales, the gravi- 
tational infall into overdense regions preferentially com- 
presses the correlation function in the a direction. There- 
fore, peculiar velocities can be statistically measured by 
analyzing the observed anisotropic pattern of £ s (a, tt) in 
both the linear and non-linear regimes. 

£ s (a, tt) is derived from the convolution of £ (r) with 
a probability distribution function of peculiar velocities 
along the line of sight, which is usually called the stream- 
ing model [25| . Even with the simplest form of a Gaussian 
probability distribution, the streaming model describes 
the suppression effect on £ s (ct, tt) on small scales. 

In the linear regime, the density fluctuations and pe- 
culiar velocity are coherently evolved through the conti- 
nuity equation, which is known as the Kaiser limit. Thus 
the known correlation function of £(r) from the linear 
perturbation theory developed by gravitational instabil- 
ity is uniquely transformed into £ s (cr, 7r) (26l - [3Cl |. 

The large scale limit of the streaming model is con- 
sistent with the Kaiser limit (22|, when both the density 
and peculiar velocity fields are treated as statistical quan- 
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FIG. 1: A£ ; *(r) for various CMB experiments priors. In the upper 
panel, we show the change in £;*( r ) for variations in u) m . In the 
lower panel, we show the change in C;*( r ) as a function of nj. The 
thick black curves are based on WMAP priors, while the thin blue 
curves are for Planck prior. On the y-axis, we focus on the range of 
scales probed by recent and planned reshift-space distortion mea- 
surements. 



tities [31|. This consistency test was developed further 
to show that, even in the Kaiser limit, the description 
of £ s (cr, tt) in linear theory can be modified due to the 
correlation between the " squashing" (in the a direction) 
and dispersion effects (in the tt direction) (32j. With 
the assumption of a Gaussian pair-wise velocity distribu- 
tion function, the dispersion effect smears into the Kaiser 
limit description of £ s (a, tt) at around the percent level 
which for our present work can be ignored. Thus we 
adopt the Kaiser limit for the description of £ s (er, tt) in 
linear regime while considering dispersion effect as a sys- 
tematic uncertainty We introduce below a new param- 
eterisation of £ s (a, tt) in terms of the cosmological pa- 
rameters and construct a method to measure the mean 
velocity dispersion v p in a model independent way. 



A. Model independent parameterisation of power 
spectra 

The discovery of cosmic acceleration has prompted 
rapid progress in theoretical cosmological research and 
prompted many authors to propose modification to the 
law of gravity beyond our solar system. For example, 
some theories based upon General Relativity can be mod- 
ified by screening (or anti-screening) the mass of gravi- 
tationally bound objects @ , while others include a non- 
trivial dark energy component (e.g. interacting dark en- 
ergy [33l . |34| , or clumping dark energy [35j ) thus break- 
ing the dynamical relations between density fluctuations 
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and peculiar velocity in the simplest dark energy mod- 
els. These theoretical ideas motivate us to express various 
power spectra of the density field in a more convenient 
way to test such theories. 

We assume a standard cosmology model for epochs ear- 
lier than the last scattering surface, and that the coherent 
evolution of structure formation from the last scattering 
surface to the present day is undetermined due to new 
physics relevant to the cosmic acceleration. Thus we di- 
vide the history of structure formation into two regimes; 
epochs before matter-radiation equality (a eq ) and a later 
epoch of coherent evolution of unknown effect on struc- 
ture formation from new physics. We can then express 
various power spectra of the density field splits into these 
two epochs, with the shape-dependent part determined 
by knowledge of our standard cosmology, and the coher- 
ent evolution part only affected by new physics. Mathe- 
matically, this is written as, 



P$$(/c, a) 
Pbb(k, a) 
p m e m (k,a) 



= £>*(fc)ff|(a), 
= D m (k)g 2 (a), 
= D m {k)g%(a) , 



(1) 



where $ denotes the curvature perturbation in the New- 
tonian gauge, 



ds z 



(1 + 2^)dt 2 + a 2 {l + 2<£)dx 2 



(2) 



and m denotes the map of the re-scaled divergence of 
peculiar velocity 9 m as m = 9 m /aH. These power 
spectra are then partitioned into a scale-dependent part 
(D$(k),D m (h)) and a scale-independent (coherent evolu- 
tion) component gb, g& m )- We define here gt = bgs m 
where b is the standard linear bias parameter between 
galaxy (or cluster) tracers and the underlying dark mat- 
ter density. 

The shape of the power spectra is determined be- 
fore the epoch of matter-radiation equality. Under the 
paradigm of inflationary theory, initial fluctuations are 
stretched outside the horizon at different epochs which 
generates the tilt in the power spectrum. The predicted 
initial tilting is parameterised as a spectral index (ns) 
which is just the shape dependence due to the initial 
condition. When the initial fluctuations reach the coher- 
ent evolution epoch after matter-radiation equality, they 
experience a scale-dependent shift from the moment they 
re-enter the horizon to the equality epoch. Gravitational 
instability is governed by the interplay between radiative 
pressure resistance and gravitational infall. The different 
duration of modes during this period results in a sec- 
ondary shape dependence on the power spectrum. This 
shape dependence is determined by the ratio between 
matter and radiation energy densities and sets the loca- 
tion of the matter-radiation equality in the time coordi- 
nate. As the radiation energy density is precisely mea- 
sured by the CMB blackbody spectrum, these secondary 
shape dependences are parameterised by the matter en- 
ergy density u> m — Q m h 2 . Both of these parameters are 
now well-determined by CMB experiments. 



D *(k) = ^^l(k)Ti(k) 



(3) 



The shape factor of the perturbed metric power spectra 
D^(k) is denned as 

2tt2_9_ 

which is a dimensionless metric power spectra at a eq , 
where A 2 a (fc) is the initial fluctuations in the comov- 
ing gauge and T$(fc) is transfer function normalized at 
T*(fc — > 0) = 1. The primordial shape A 2 g (k) depends 
on ns, as A 2 (k) = A 2 s (k/k p ) ns ~ 1 , where A 2 S is the am- 
plitude of the initial comoving fluctuations at the pivot 
scale, k p = 0.002Mpc _1 . The intermediate shape factor 
T$(fc) depends on oj m . The shape factor for matter fluc- 
tuations and peculiar velocities D m (k) are given by the 
conversion from D<&(k) of, 

D m(k) = tj^D^k), (4) 

where = 1/2997 Mpc" 1 . 

Unlike the shape part, the coherent evolution compo- 
nent, gb and ge m are not generally parameterized 
by known standard cosmological parameters. We thus 
normalize these growth factors at a eq such that, 

g^{a eq ) = 1, 
gs m { a eq) = a e g3*(a eq ), 

dgS m ( a eq) 



9e m {aeq) = — 



dlni 



(5) 



Instead of determining growth factors using cosmologi- 
cal parameters, we measure these directly in a model- 
independent way at the given redshift without referenc- 
ing to any specific cosmic acceleration model and with 
the minimal assumption of coherent evolution of modes 
after a eq . Considering the uncertainty in the determi- 
nation of A 2 S from the CMB anisotropy, which is degen- 
erate with the optical depth of re-ionization, we com- 
bine both A 2 S and gx (where X denotes each compo- 
nent of <E>, b and m ) with proper scaling for conve- 
nience as g* x = gxAs/Ag. Throughout this paper, we 
use A* s 2 = 2.41 x 10~ 9 for the mean A s value from the 
WMAP5 results. Our result on measuring the bulk flow 
motion is independent of our choice of an arbitrary con- 
stant A* s 2 . 



B. Correlation function in the configuration space 

In the linear regime of the standard gravitational in- 
stability theory, the Kaiser effect (the observed squeezing 
of £ s (<t, 7r) due to coherent infall around large-scale struc- 
tures) can be written in configuration space as, 

£ s faOfa) = U 2 + ^.9e m + J.9e 2 m )eo(0^o(/i) 
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(6) 
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FIG. 2: The £ s (cr, w) correlation function. We plot for three con- 
tours of § s (ct, 7r) = 1,0.1,0.01 (from the inner to outer contour). 
In the upper panel, the solid curves are for a ACDM cosmology, 
while the dash and dotted curves are for models with g b = 1.5 and 
2 respectively. In the lower panel, the solid curves are for a ACDM 
cosmology, while the dash and dotted curves are for models with 
3e m = 1.5 and 2 respectively. 



where Pi (/i) is the Legendre polynomial and the spherical 
harmonic moment £* (r) is given by, 



k 2 dk 



(7) 



where ji is a spherical Bcssel function and * denotes scal- 
ing of the shape factor with A* s 2 . 

As discussed above, we ignore the effect on £ s (er, it) of 
the small-scale velocity dispersions within a single dark 
matter halo [32| as the effect is only a few percent, and 
split Eqn. |B]into a shape-dependent part which 
is determined by the cosmological parameters (n$ and 
w m ), and a coherent evolution component, which is pa- 
rameterised by g^, <?g m i at the targeted redshift z.;. The 
shape part is therefore almost completely determined by 
CMB priors, while the coherent evolution of structure 
formation can be determined from fitting £ s (<x, it) , in 
redshift -space, as a function of redshift. 

In Figure [TJ we present the effect of CMB priors on 
the value of £*(»")■ I R the top panel of Fig. [TJ we provide 
the expected variation in £f(r) from varying uj m . We see 
that varying a eq causes greater tilting in the shape of 
££(r), since larger scale modes can come into the hori- 
zon earlier. In addition to this contribution, the overall 
amplitude of £* (r) depends on uj m by a weighted trans- 
formation between D$(k) and D m (k). Considering the 
marginalisation over CMB priors, we expect a discrep- 
ancy of ~ 5% with WMAP5 measurements, and just a 



few percent effect with the projected Planck priors. 

In the bottom panel of Fig. [1] the dependence of £ ; * (r) 
on ns is given for both WMAP5 and Planck priors. The 
overall shifting can be re-scaled by adjusting the pivot 
point to the effective median scale of the survey. With 
the measured WMAP5 prior of Ans = 0.015 [36|], we ex- 
pect variations of a few percent on the sha pe, while for an 
estimated Planck prior of Ans — 0.0071 [33|, we expect 
£;*(r) to be nearly invariant to ns- The shape of £f (r) is 
affected maximally during the intermediate epoch, from 
horizon crossing to the matter-radiation epoch. The de- 
cay rate of the inhomogeneities differs by the ratio be- 
tween matter and radiation energy densities. 

Once CMB constraints are placed on the shape part 
of Eqn. 6, the coherent history of structure formation is 
obtained from the anisotropic moment of £ s (<7, 7r). Even 
though both and ffe m i weight the evolution sector 
simultaneously, their contribution to £ s (a, it) are differ- 
ent, which enables us to discriminate g^ i from <^ -• In 
the monopole moment, g£ li is the dominant component 
since g£ >i > gQ m , unless their is an excessive bulk flow. 
Thus the variation of <j£ i generates a near isotropic am- 
plification as illustrated in the top panel of Figure [2] In 
the quadrupolc moment, the cross-correlation between 
5 m and 6 m is leading order. The reversed sign of the 
quadrupole moment results in the squashing effect, and 
it is sensitive to the variation of c/q i as the cross- 
correlation is the leading order. In the bottom panel 
of Figure [2] the variation of <j@ i mainly contributes to 
the anisotropic moment. It is this signal which allows 
both g£ )i and <?q i to be probed separately using the 
anisotropic structure of £ s (<7, it). The contribution from 
the term having peculiar velocity autocorrelation is not 
significant if excessive bulk flows are excluded. 



C. Implication for cosmology from measuring g1 i 

and 3e m * 

A measurement of g^ m i is equivalent to the quantity 
y (J mass JJJJ anc j th ere f ore ; an excellent test of dark en- 
ergy models (where / is the logarithmic derivative of 
the linear growth rate and cr™ ass is the root-mean-squarc 
mass fluctuation in spheres with radius 8/i _1 Mpc). While 
cosmological test of g^ . are free from bias, which is 
notoriously difficult to measure accurately in a model- 
independent way, the reported value of g^ m i does depend 
on the normalization which is also poorly constrained 
(i.e., primordial amplitude) or model-dependent (i.e., 

_mass \ 
°8 )■ 

Thus, we introduce a more convenient paramctcrisa- 
tion of peculiar velocity which is independent of these 
normalization issues. The measured g^ i (in the redshift 
bin i) which can be translated into the one-dimensional 
(1-D) velocity dispersion in that redshift bin (u*) by, 
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FIG. 3: The 2-D two— point correlation function (£s(c, 7r)) for the 
SDSS DR3 C4 cluster survey with a median redshift of 2=0.1. The 
contours have been slightly smoothed. 



In this formula, there is a degeneracy between g? t and 
3q i which cannot be solely broken by fitting £ s (cr, 7r); 
instead we simultaneously fit for vt and b l from the data 
and then marginalize over the bias to obtain v p (indepen- 
dent of b) in that redshift bin. Therefore, if our statistical 
determination of the history of v p can be combined with 
an independent measurement of bias, then v l p can be de- 
termined precisely. The scaled parameter c/g, m i depends 
on all shift factors; the primordial amplitude or the en- 
hancement of D m due to varying u m , as well as later 



evolution. But the estimation of v p from <?g 



time O r 

is independent of the uncertainty in the overall shifting. 

If the evolution of g^ is measured, it can be used to 
reconstruct other coherent growth factors. The coherent 
growth factor of $ can be given using the Euler equation, 



2 aH 



9<t> = 



3 H*UJr, 



9e„ 



d 9*e 



d\n a 



(9) 



where no anisotropy condition is used. If the Poisson 
equation is validated then the re-constructed g^, can be 
used to derive using the relation g* s = ag^. Then 
this estimated matter fluctuation evolution can be used 
to determine bias from the measured as b = gli/gg m - 



III. REDSHIFT— SPACE DISTORTIONS FROM 
CLUSTERS OF GALAXIES 

As a demonstration of the parameterization discussed 
above, we present here a measurement of the redshift - 
space 2-D two-point correlation function (£ s (<7, 7r)) for 
clusters of galaxies selected from the SDSS. We use an 
updated version of the C4 cluster catalogue [23| based 



on Data Release 3 (DR3; [3j|) of the SDSS. Briefly, the 
C4 catalogue identifies clusters in a seven-dimensional 
galaxy position and colour space (righ ascension, declina- 
tion, redshift, u — g,g — r,r — i,i — z) using the SDSS Main 
Galaxy spectroscopic sample. This method greatly re- 
duces the twin problems of projection effects and redshift 
space distortions in identifying physically-bound galaxy 
groups. This catalogue is composed of ~ 2000 clusters in 
the redshift range 0.02 < z < 0.15. 

The estimation of the correlation function relies cru- 
cially on our ability to compare the clustering of the data 
to that of a random field. Thus any artificial structures 
in the data must be considered when constructing the 
random catalogue. These problems include incomplete- 
ness, such as the angular mask (e.g. survey boundaries, 
bright stars and dust extinction in our own galaxy), and 
the radial distribution where at large distances, the mean 
space density decreases as we approach the magnitude 
limit of the survey. We have constructed random sam- 
ples which takes these issues into account, i.e., the angu- 
lar positions are randomly sampled from a sphere to lie 
within the DR3 mask, while the redshifts are obtained 
from a smooth spline fit to the real C4 redshift distri- 
bution (which removes true large scale structures). The 
random samples are then made to be 50 times denser 
than the real data to avoid Poisson noise. 

In Figure [31 we show our estimation of the £ s (a, n) 
binned into with 6 configuration-space bins up to 60 Mpc 
(one bin per lOMpc). Separations of less than 10 Mpc 
are removed to reduce the FoG effect. Error on £ s (er, 7r) 
were derived using the jackknife method [3{|, which in- 
volves dividing the survey into N sub-sections with equal 
area (and thus volume) and then computing the mean 
and variance of £ s (c, ir) from these N measurements of 
the correlation function with the i th region removed each 
time (where i = 1...N). 

In our analysis, we divided the whole C4 area into TV = 
30 sub-subsections and determine the variance from 14011 , 
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4(n) = %^ £[&(n)-ec 
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n)\ 2 , 



jack 



(10) 



fe=i 



where Nj ac k is the number of jackknife samples used and 
n represents a single bin in the a — tt configuration space. 
Then we compute 



Njack 



jack 



(ii) 



fe=i 



and the normalised covariance matrix is estimated 
from 39], 



1 v jack 



k=N jack 
k=l 



where, 



.fe gfc(n) -£(rj) 

"tin) ' 



(12) 
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FIG. 4: An ordered list of the Eigenvalues for our cluster covari- 
ance matrix. 



Before we invert in Eqn.[12] we note that the values 
of Cij are estimated to limited resolution, 




(14) 



and therefore, if Njk is small, or there are degeneracies 
within Cij, the inversion will be affected. This problem 
can be eliminated by performing a Single Value Decom- 
position (SVD) of the matrix, 



Cij = UlD kl V ljt 



(15) 



where U and V are orthogonal matrices that span the 
range and the null space of Cy and D^i = A 2 <5fc/, a diag- 
onal matrix with singular values along the diagonal. In 
doing the SVD, we select the dominant modes t o con- 
tribute to the x 2 by requiring that A 2 > y/2/Njk. 

In Figure HI we rank the eigenvalues (Aj) for the in- 
creasing cigenmodes and see a "kink" in the distribu- 
tion which we interpret as indicating a transition in the 
signal-to-noisc of the cigenmodes, i.e., only the first ten 
modes contain most of the signal, while higher-ordered 
modes are dominated by noise. Wc therefore remove 
cigenmodes beyond this kink (with A» < 0.01) where the 
eigenvalues start to flatten out. 



A. Statistical determination of large scale flow 

As discussed in the Introduction, there is recent evi- 
dence for excessive bulk flow motions compared to the 
WMAP5-normalised ACDM model [3 and therefore, it 
is important to confirm these results as it may indicate 



evidence for an alternative explanation for the observed 
cosmic acceleration such as modified gravity. In this pa- 
per, we provide a first demonstration of our new param- 
eterization using clusters of galaxies from the SDSS. In 
detail, we attempt to model the "squashing" of the 2- 
D correlation function of the C4 cluster sample seen in 
Figure 5 using the formalism presented herein. We do 
however caution the reader that we expect the limited 
size of the DR3 sample to leads to large statistical er- 
rors, due to a significant shot-noise contribution because 
of their low number density. However, future cluster and 
galaxy samples (e.g., LRGs) should provide stronger con- 
straints and provide a more robust test of these high bulk 
flow measurements in the literature. 

In Figure [5J wc provide the best fit parameters b and 
v p for the C4 correlation function presented in Figure 
3 and there is as expected a clear anti-correlation be- 
tween these two parameters because the anisotropic am- 
plitude is generated by cross-correlations in the density 
and peculiar velocity fields. The best fit value from Fig- 
ure 5 is v p = 270 +433 km/s (at the la level marginalised 
with other parameters including b) and is consistent with 
v p = 0. Wc do not quote the negative bound of the error 
on v p as it is below zero and thus has no physical mean- 
ing. Instead, we quote the upper bound on v p and note 
that our result is consistent with zero. Our measurement 
of v p is close to the predicted value of 203 km/s for a 
WMAP5-normalised ACDM model. 

We propose above that v p is a complementary param- 
eter for reporting such peculiar velocity measurements. 
The parameter g@, which is equivalent to /os, is not de- 
termined precisely without the prior information of As- 
But when we report our measurement with v p , there is no 
uncertainty due to other cosmological parameters which 
are not determined accurately, as it is equivalent to g^ 
determined statistically from redshift space distortion. 
The observed value v p at a given redshift is not only in- 
dependent of bias but also independent of normalisation. 



B. Reconstruction of matter density field from v v 

We convert v p measurement into go using A s from 
WMAP5 (ge : coherent growth factor of peculiar veloc- 
ity, and it is equivalent to fag in other parameteriza- 
tions). With the evolution of <7e known, dynamics of per- 
turbations are reconstructed to provide the history of \t 
through the Euler equation. In most theoretical models, 
the time variation of v p is minimal at these low redshifts 
discussed here for the C4 sample (z ~ 0.1), which allows 
us to ignore the time-derivative part in Eq. [51 Therefore, 
it is straightforward to transform the coherent evolution 
of into the coherent evolution of ty. If we assume no 
anisotropic stress, then it is easy to convert to the coher- 
ent evolution of $, g$. 

We are able to determine matter density fluctuations 
through the Poisson equation. We calculate the coher- 
ent growth of 6 m , gs = 0.7, which is related to g$, as 
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FIG. 5: The 2-D contours between b and v p with the DR3 cluster 
sample. 



gs = ag®, if no modified Poisson equation is assumed. 
Finally, the estimated gs can be used to derive bias using 
measured g£. Through fitting to the redshift distortion 
effect, we extract both g£ and g^. The density fluctu- 
ation evolution gg is estimated only from <?q, and the 
other measurement g£ is not yet used. The combination 
of the estimated gs and the measured g@ provides bias 
from b = gl/gs = 2.9 ±0.8, which is fully consistent with 
our expectations for such massive clusters of galaxies in 
the C4 sample. 

IV. DISCUSSION 

We outline in this paper a new theoretical model 
for £ s (7r, cr), the 2-D two-point correlation function in 
configuration-space, which allows us to constrain the 
bulk flow motion of matter on large scales. We also 
propose that the 1-D linear velocity dispersion (v p ) is a 
interesting quantity to report when measuring rcdshift- 
space distortions. We demonstrate this method using C4 
clusters from the SDSS and find a value for v p that is 
consistent with a WMAP5-normalised ACDM cosmol- 
ogy (within our large statistical errors). Our observed 
value for these bulk flows is marginally inconsistent with 
other recent observations in the literature, which find an 
excess flow compared to a WMAP5-normalised ACDM 
model fl8l - l2"oT | . We do not discuss this further as we plan 



to revisit these measurements using larger datasets and 
different tracers of the density field. 

As discussed in Section II- C, our measurement of v p 
is correlated with gt,, which is the combination of b and 
6, n , since the observed anisotropic shape of the 2-D cor- 
relation function in redshift-space is generated by a cross 
correlation between the density field and peculiar veloc- 
ities. However, one of the important implications of our 
method is that we can measure v p without knowing how 
to decompose <?{,, and thus without the uncertainty of 
determining b. 

There are however some caveats to our analysis. We 
do not analyse our data in Fourier space, but in config- 
uration space. Small scales have been removed from our 
data analysis (< lOMpc), to ensure that the FoG effect 
will not contaminate our results. Our methodology is in- 
sensitive to a possible shape dependence at large scales 
for any exotic reason; scale dependent later time growth 
(e.g. f(R) gravity models [4lj). or scale-dependent bias 
at large scales. In follow-up studies, we will measure 
the redshift-space distortions in Fourier space to test the 
effect of small-scales on our results. In addition, the for- 
mulation to derive £ s used in this paper can be slightly 
biased due to the dispersion effect studied in [32[ . This ef- 
fect is not parameterised properly here, but the reported 
level of uncertainty is approximately 5% which is much 
smaller than the statistical errors on our present mea- 
surements. Therefore, we dismiss this shift here but it is 
worth revisiting this issue in the future to know how to 
incorporate this effect in a new parameterisation. 

Finally, Song and Percival [13[ recently proposed a 
method to re-construct the structure formation observ- 
ables from O measurements. Although it is not yet es- 
timated precisely in that narrow range of measured val- 
ues, we apply their methodology in practice. From the 
observed coherent evolution of O at z = 0.1 from the 
DR3 C4 clusters, we re-construct $ and 8 m . We then 
find that bias can be derived from the estimated S m and 
the measured g^. Here, for the first time, we estimate 
bias from peculiar velocity measurements only. The esti- 
mated values are reasonable at b = gl/gs = 2.9 ± 0.8. It 
is not precise measurement yet, as the time variation is 
ignored, but we will revisit this in a following paper. 
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